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Abstract 

We show that a right-handed stop in the 200-400 GeV mass range, 
together with a nearly degenerate neutralino and, possibly, a gluino 
below 1.5 TeV, follows from reasonable assumptions, is consistent 
with present data, and offers interesting discovery prospects at the 
LHC. Triggering on an extra jet produced in association with stops 
allows the experimental search for stops even when their mass dif- 
ference with neutralinos is very small and the decay products are 
too soft for direct observation. Using a razor analysis, we are able to 
set stop bounds that are stronger than those published by ATLAS 
and CMS. 



1 Introduction 



Supersymmetry has been significantly cornered by LHC searches. The discovery of the 
Higgs boson at 126 GeV [1, 2] and the direct limits on sparticles rule out most of the natural 
implementations of low-energy supersymmetry, at least in their simplest versions [3]. Pockets 
of parameter space still survive, but their exploration requires the 14-TeV phase of the LHC. 
At this stage, it is appropriate to examine the available experimental data and look for hints 
that can guide us towards special regions where supersymmetry may still hide. 

In this paper, we point out that there is a window of supersymmetric parameters that 
(i) are well consistent with all collider data and flavor constraints, (ii) naturally emerge 
from RG evolution of simple UV completions, (in) predict the correct thermal abundance 
for dark matter (DM), and (iv) give observable signals at LHC14. In this special window, 
the supersymmetric mass spectrum has the following properties: 

• The lightest stop is mostly right-handed and its mass is in the range mi = 200- 
400 GeV. 

• The heavy, mostly left-handed, stop has a much larger mass (in the 1-2 TeV range), 
but it is correlated with the light stop in such a way that their geometric average is 
m s = {m- H m h fl 2 « 500-600 GeV. 

• The stop trilinear term is large, such that 1 A\ « 6m|. 

• The gluino mass is below about 1.5 TeV. 

• The lightest neutralino has a mass slightly smaller than the lightest stop, by an amount 
of about 30-40 GeV. 

In section 2 we will give several arguments that lead to the mass spectrum described 
above. None of them is sufficiently compelling to select conclusively the sparticle masses 
but these arguments, taken together, give circumstantial evidence in favor of our choice of 
parameters. Our conclusion is based on the following considerations: 

• We choose the values of the stop parameters that minimize the average stop mass, while 
leading to a Higgs mass of about 126 GeV. The resulting stop mass spectrum, although 
not strictly natural, has the advantage of reducing the amount of unnaturalness forced 
upon supersymmetry by present LHC data. 

1 This configuration is known as the "maximal mixing" case, although it does not necessarily imply a 
large mixing between the two stop mass eigenstates, as discussed in sect. 2. 
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• Stops affect the rates for gg — >■ h and h — > 77. Experimental data are at present 
not sufficiently constraining, but will soon play an important role in selecting stop 
parameters. In particular, we find that stop contributions to gg — > h and h — > 77 
exactly cancel for m~ t2 m 6mj , under the conditions preferred by the Higgs mass value. 

• Flavor constraints rule out a light left-handed stop, but are consistent with a light 
right-handed stop. We also show that, if the CKM matrix is the only source of flavor 
violation and the higgsino is relatively light, supersymmetric contributions can improve 
the agreement with the measurements of e#, while being compatible with B — > X a j. 
This is achieved for a large mass splitting between left-handed and right-handed stops. 

• A large splitting between left-handed and right-handed stop masses naturally emerge 
from RG evolution, as long as the gluino is not too heavy. Moreover, once we require 
a large stop mixing parameter at low energy, we find an upper bound on the gluino 
mass. 

• A light stop can be very helpful to obtain the right dark-matter relic density, which 
is typically too large for B-ino LSP or too small for higgsino or W-ino LSP in generic 
supersymmetric models. The process of coannihilation selects the preferred value of 
the mass difference between stop and neutralino. 

After we have determined the favorable region for sparticle masses, we study in sect. 3 
the experimental strategies for discovery at the LHC The phenomenology of supersymmetric 
frameworks with light stop has been discussed at length in the recent literature (see e.g. ref. [4, 
5, 6, 7, 8]). The search of light stops has also been the focus of recent dedicated experimental 
analyses by both ATLAS and CMS [9]. However, the peculiarity of our scenario is the 
near mass-degeneracy between stop and neutralino, which makes experimental identification 
especially arduous (see ref. [10, 11, 12] for previous attempts to address this problem). 

On the one hand, we show that the decay t — > Nb£i>e, which has been neglected in 
the present exclusive experimental searches, can dominate over the more traditional decay 
t — > cN, especially if the mass difference between t and N is not too small. The four-body 
decay process has the advantage of producing observable leptons in the final state, leading 
to a possibly higher signal/background ratio in exclusive searches. On the other hand, by 
an explicit simulation of this decay channel and the analysis of presently available data, we 
show that the inclusive searches, and particularly the CMS razor analysis, already provides 
significant constraints on this framework. 
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Figure 1: The Higgs mass in low-energy super- 
symmetry for large tan/3 ~ 20. The shaded re- 
gion in the (Xt,ms) plane corresponds to the 
observed value of m^. Higher-order corrections 
and the uncertainty in the top mass amount to 
an error of a few GeV in m,h . 

2 The light-stop window 



Figure 2: The white region is the range in the 
(m£ , mj ) plane allowed by the nth constraint, 
while shaded regions are excluded. The full, 
dashed, and dotted lines correspond to fixed val- 
ues of At, satisfying the constraint with 
\X t \ > V6 (blue) or \X t \ < y/d (black). 



2.1 Constraints from the Higgs mass and decay rates 



The leading part of the supersymmetric prediction for the mass of the lightest Higgs boson 
is 



m. 



m 2 z cos 2 (3 + 
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+ 
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where X t = (At + \i cot (3)/ms, m 2 s = mf is the average stop mass, y t = m t /v is the 
top-quark Yukawa coupling, and v ~ 174 GeV is the Higgs vev. In fig. 1 we show the region 
of the (X t ,ms) plane compatible with the observed Higgs mass (for tan/3 3> 1), including 
also the leading two-loop corrections to the Higgs mass not shown in eq. (1). The lightest 
average stop mass that can lead to the observed Higgs mass is obtained for 



m s w 500 GeV and X] 
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(2) 



We focus on such configuration, the so-called "maximal mixing" case, since it reduces the 
fine-tuning in electroweak symmetry breaking and can lead to observable signals. 
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Further constraints on X t and the stop masses can be obtained by examining the correc- 
tions to the h — > 77 and h <H- gg rates: 

r^TT =(1 + A) , r(A->^) =(1 _ oWi (3) 
r(n ^ s-^jsM r(/i -»■ 77)sm 

where, in the limit in which we decouple the pseudoscalar Higgs, we find 

Present data (fitted in the context of the SM plus light stops) give [13] 

A t = -0.04 ±0.11 (5) 

and do not yet imply a significant constraint, as it is clear from fig. 2 where we plot iso-curves 
of A t after imposing the requirement. The situation will improve in the future. Note 

that no deviations from the SM (A t m 0) are obtained for m$ ~ 6777.^ if we insist on having 

2 
t 

A few comments are in order: 



An independent indication of a large splitting between m t - 2 and mi can be obtained if 
we assume that A t is not significantly larger than the trace of the stop mass matrix. 
Assuming A\ < a(m~ + m| ), then (for large tan/3) X t 2 is bounded by 



ml + ml r<<1 a 



-1^9 ti ' to r<J f- i " "Hi /n\ 

X: < a—^ ^ - - , r = . (6) 

m h m- t2 r m h 

Vacuum stability arguments imply a < 3 (assuming m 2 Hu <C ?« 2 2 ), but this does not 
allow us to deduce a significant constraint on r. However, if a £ 1 (as naturally 
expected from RG arguments, see next section) then we are forced to assume small 
values of r in order to reach X 2 w 6. 

Despite the large value of X t , the mixing of the two stop eigenstates is suppressed in 
the limit r < 1: 

1 . / 2m t m s X t \ r «i rliWit 

t = - arcsm — = =- ~ . (7) 

2 \ m i 2 ~ m i j m s 

So, in this limit, we can approximately identify the two mass eigenstates with the 
electroweak eigenstates. As we will show in the next section, it is natural to identify 
the lightest state with an almost right-handed stop. Note also that for r 1 the 
lightest stop mass is significantly lighter than the average stop mass in eq. (2): r 1/6 
corresponds to m^ ~ 200 GeV. 
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Figure 3: Renormalization group evolution from 
the unification scale to the weak scale of gaug- 
ino masses M\, M 2 , M3 (green curves), of 
the stop mass parameters and mj^ (full 
and dashed blue curves, respectively), ytA t (red 
dashed curve), rriH u (black curve), in a configu- 
ration leading to m^ R <C m i L a t the weak scale. 
All masses are in GeV units and we assumed 
the MSSM. 



Figure 4: Gluino and light-stop masses result- 
ing from a scan of the parameter space assum- 
ing universal scalar and gaugino masses, and 
the condition \At\ < 3mo, at the GUT scale. 
All points satisfy the rrih ~ 126 GeV con- 
straint and are colored according to the value of 
mj /mj , as indicated on the right-handed axis. 
For illustrative purposes lines corresponding to 
M^/mr = 1,2,3,4 are also shown. 



2.2 Constraints from the RG evolution 



A numerically large splitting between m t - and rrii R naturally arises from the evolution under 
renormalization-group equations (RGE), provided scalar masses are significantly larger than 
gaugino masses at the high scale. This can be understood by looking at the one-loop RGE 
for third generation squark masses and m# u . Neglecting off-diagonal flavor-mixing terms we 
have 

8- 2 ^ = V 2 t Y t - f 9 M - 3g 2 Ml - l^glMl , (8) 

= ^r t -f 9 M-f 59 M, (9) 

Svr 2 -^ = (10) 
a log fi 5 

where ft is the renormalization scale, and 

Yt = ml+ml R + m 2 Hu + Al . (11) 
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The RG evolution of the stop masses depends mainly on two effects: the QCD term (g|M|) 
and the Yukawa term (y 2 Y t ). If we take M 3 , m~ th and A t to be comparable and in the range 
1-2 TeV at the weak scale (in order to fulfill the m^ constraints and the experimental bounds 
on the gluino mass), we find that: i) QCD and Yukawa terms compensate to a large extent 
in the running of mi L ; ii) the Yukawa term is dominant during most of the running of mf , 
leading to mi R <C mi L at the weak scale starting from the initial condition m^ R = m^ at 
some high scale; iii) the Yukawa term is always dominant in the running of m? H , which 
naturally becomes negative at the weak scale. 

An illustrative spectrum is shown in fig. 3, where we required rn 2 Hu = — m|/2, M 3 = 
1.3 TeV, and m t - < 300 GeV at the weak scale, and adjusted A t in order to achieve the 
condition mi R = mi L at 2 x 10 16 GeV. The corresponding weak-scale configuration is consis- 
tent with all the existing experimental bounds, with the condition rrih = 126 GeV (assuming 
tan/3 > 5), and with a light stop below 300 GeV. A few percent tuning in the initial values 
of m 2 H and the higgsino mass ix is necessary in order to achieve the correct pattern of elec- 
troweak symmetry breaking, but this is unavoidable in the minimal supersymmetric model 
with rrih = 126 GeV. The soft-breaking terms needed to reach this low-energy configuration 
require an initial splitting ra 2 quar ks/ m gauginos ~ 10 ^ the high-energy scale. All three genera- 
tions of squarks can be degenerate at the high scale, since the separation of the right-handed 
stop is fully driven by the dynamics of the low-energy degrees of freedom. 

The dynamical separation between m^ R and m; L , together with the generation of a large 
X t , from high-scale RG running naturally occurs only in a limited range of gluino masses. 
This can be understood by inspecting the expressions of A t , mi R and raj at the weak scale 
in models with a universal scalar mass mo and trilinear coupling Aq at the GUT scale, 

A t « 0.3A + O.8M3 , (12) 
m\ R w 0.5M 3 2 - 0.07A 2 - O.IOA0M3 + 0.3m 2 , (13) 
m\ h w 0.7M 3 2 - 0.03A 2 - 0.05A M 3 + 0.7m 2 , (14) 

where M 3 is the gluino mass at the weak scale. From these equations it is clear that if 
M3 <C mo, |Aq| a large splitting among the two stop masses and a large X t are obtained only 
for unnaturally large values of |y4 |/m . Similarly, maximal mixing and large splitting cannot 
be obtained if M 3 3> mo, |-4o|- The upper bound on M 3 , which is particularly important 
for the LHC searches, is quantified in fig. 4, where we show the points satisfying the rrih 
constraint in the m^-M 3 plane. As can be seen, the gluino must satisfy the approximate 
upper bound M 3 ^ 4m^, that implies M 3 ^ 1.6 TeV for the range of m^ (m^ ^ 400 GeV) 
corresponding to a large rn^/m^ ratio. 
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As anticipated, in this framework m 2 H naturally becomes negative at the weak scale 
and the \x term must be properly adjusted to reproduce the correct value of mz- Assuming 
universal scalar masses at the high scale, m% runs very negative at the week scale, implying 

> mi R) or heavy higgsinos. Alternatively, we can consider a scenario as in fig. 3, with 
non-universal boundary conditions, where \y\ < m^ R and thus higgsinos are lighter than the 
right-handed stop. The two cases lead to a rather different phenomenology for flavor, dark 
matter, and LHC searches. 



2.3 Constraints from flavor physics 

If some of the gauginos or higgsinos are not too heavy, a light stop can have a significant 
impact on low-energy flavor-physics observables. On general grounds, even if gauginos and 
higgsinos are in the several- TeV domain, sizable misalignments in flavor space between quark 
and squark mass matrices are excluded. Therefore we assume that the light stop is mostly 
right-handed and aligned in flavor space with the top quark. The remaining flavor violation 
is described by the usual CKM angles in charged currents. An interesting and largely model- 
independent correlation (controlled only by the size of and the stop mass parameters) 
emerges between BR(B — > X s ^) and ex- 

In the limit in which we retain only the effect of higgsinos and flavor-aligned right-handed 
stop, the deviations from the SM in these two observables are described by 

BR(B^X sl ) 1 _ 2 x ACr _M ACai (15) 



BR(B X sl ) 



SM 




(16) 



2 \ ™2 



AC 7 , 8 = sin 9 t tan fl^F% MU + Ir^s 4" > ( 17 ) 




where 

mi ' \mj I mi 

and the normalization of the various loop functions (see Ref. [15, 14] for the explicit expres- 
sions) is 

F 7 LR (1) = -- , F 8 LR (1) = -— , F 7 RR (1) = — , F RR (1) = — , F 2 (l) = — . (18) 

7 V ) g , 8 V J 12 » 7 V ) 2g8 ) 8 V J gg » 12 V 7 

For B — > X s 7 we have expanded the result to first order in the stop-mixing angle. Note 
that, even for | sin# t | 1, the first term in eq. (17) can be sizable and can dominate over 
the second one, both because of the large value of the loop function Fj R and because the mu 
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Figure 5: Correlation between BR(I? — > X s ^f) and ex- The two ellipses denote the 68% and 90% CL 
experimental range. All points reproduce the observed Higgs mass. The two black curves are obtained 
varying between 200 GeV and 400 GeV (from left to right) for ms = 500GeV ; fi = 250 GeV, 
and tan/3 = 20 (dashed curve) or tan/3 = 10 (full curve). The points are obtained varying the 
parameters in the range \i = [150 — 400] GeV and m^ R = [200 — 400] GeV, with ms < 700 GeV and 
tan/3 = 20 (red) or tan /3 = 10 (blue). 

constraint favors tan (5 ^> 1. As a result, the experimental constraint on BR(i? — > X s '-f) puts 
a very stringent bound on the maximal value of \6t\ for higgsino masses of 0(mi R ), providing 
a further argument in favor of a sizable hierarchy between the two stop mass eigenstates [see 
eq. (7)]. The sign of the correction can be positive or negative, depending on the relative 
sign of /i and A t . The experimental data favor a constructing interference with the SM 
amplitude: BR(B X s -f) exp /BR(B X s7 ) SM = 1-09 ± 0.11. 2 

In the case of ex, the correction is always positive and, in first approximation, is inde- 
pendent from the mixing angle. As a result, the present experimental constraint e^ F '/e^ 1 = 
1.14 ± 0.10 [19] can be better satisfied if /i is not too heavy. 

The correlation between the two observables is shown in fig. 5, where we restrict the 
attention to the value of sgn(/i/l t ) favored by BR(£> — > X s j). As can be seen, after imposing 
the rrih constraint and requiring <^ 400 GeV, present data favor the configuration with 
mi R mi that maximizes the correction to €k and minimizes the impact in BK(B — > X^). 

2 This ratio is evaluated using the SM estimate from Ref. [16], and a naive average of the HFAG result 
and the latest Babar result [18] on BK(B — > X s j) cxp . 
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Figure 6: Points in the supersymmetric parameter space that lead to the correct DM abundance. 
The shaded band is the analytic approximation of eq. (22). 

2.4 Constraints from dark matter 

A light stop offers the opportunity of curing the excessive relic abundance of B-ino LSP, 
generally encountered in supersymmetric models. Indeed, the DM cosmological abundance 
can be reproduced with a B-ino thermal relic that co-annihilates with stops if 

m h = M DM + AM with AM w (30 - 40) GeV. (19) 

The relatively small mass difference arises imposing that the average annihilation cross sec- 
tion equals 

av cosmo = (2.3 ± 0.1) x 10" 26 cmV 1 (20) 

at the freeze-out temperature Tf pa Mdm/25. The dominant annihilation process is s-wave 
stop annihilation into gluons (annihilation into quarks is p-wave suppressed): 



7 4 

ct(M* gg)v = 9 ^ r . (21) 



fi 

Averaging over the components of the DM system (t and t* have 3 colors each, and the 
neutralino has 2 polarisations) we get 

,3 



av. 



-^l^-aiUl gg)v. (22) 



The band where the DM abundance reproduces the cosmological values within 3 standard 
deviations, according to eqs. (20) and (22), is shown in fig. lb. The dots represent points 
in supersymmetric parameter space which give a correct DM abundance, derived from full 
numerical simulation of the relevant Boltzmann equations (from [21]). 
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3 Experimental signals 
3.1 Stop decay rates 

Beside gluino production, the most characteristic signal comes from the (mostly right- 
handed) light stop. Dark matter considerations motivate the searches for a stop that is near 
degenerate with the neutralino LSP, with a mass difference AM = — Mdm ~ 30 GeV. 
In this configuration, the stop is usually assumed to decay according to t\ — > cN. Here we 
point out that four-body stop decays (not suppressed by flavor-changing neutral currents) 
can easily become competitive with the two-body flavor-violating decay. In the limit of small 
AM, the relevant stop decay widths are 

= 2 g Wfa. W = cm _, le, V / AM y /400GeV\ 
v ; V-Km k V 10 V V30GeV/ \ m k ) 

x 3 gHan 2 6 W AM S / AM \ 8 /400GeV\ , , 



as well as 



r(fi -> iV&ud) ss r(ti -> iV6cs) sa 3r(fi ->■ Nbtui) i = e, r . (25) 



For the decay t\ — > cN, the parameter 6 tc is the effective stop-scharm mixing angle. In 
general, 9 tc is a free unknown parameter, since it depends on the flavor structure of the soft 
terms. Assuming that it vanishes at some high scale Auv, a non zero value is generated 
by RGE effects due to the SM Yukawa couplings (even in absence of other sources of flavor 
violation) [22]. In our scenario, where t\ w Ir, the leading effect comes from an induced 
tR—c.L mixing, which can be estimated as 

.mfv ._, yty 2 b VcbV t l vA A uv = 5 / 2 TeV \ / logA uv /m \ Aan/3\ 2 

tc 16tt 2 m 2 g m Vm/V 30 / V 10 / ' lJ 

where we have omitted 0(1) loop functions depending on mass ratios of heavy squarks and 
charginos, whose average mass is denoted by rh. 

The t\ — > Nb£ + vi decay receives contributions suppressed by heavy sparticles or mediated 
only by virtual SM particles. We here focus on the latter contribution, which is dominant 
in our case. This leads to eq. (24), whose derivation is given in the appendix, together with 
the matrix element relevant for implementation in Monte Carlo codes. 

The two decay channels can dominate in different regions of the parameter space and 
become roughly comparable for 9 tc ~ 10~ 5 and stop-neutralino mass differences motivated by 
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DM considerations. However, the large model dependence of 6t c prevents us from making any 
firm conclusion. The four-body decay has a much steeper dependence on AM and becomes 
less relevant for very small AM. The decay t\ — > Nbi + v^, not previously considered in the 
literature, is interesting from the experimental point of view since it leads to an additional 
soft lepton. 

Since the t% — > cN decay only produces an unobservable soft jet and its signatures have 
been previously studied, in the following we focus our attention mainly on the four-body 
decay channels in eqs. (24)-(25), which we assume to be the dominant decay modes. As 
we discuss below, the bounds we are able to derive at present from existing LHC searches 
are largely independent from this assumption; however, the presence of a lepton in the final 
state could possibly lead to stringent bounds with future optimized searches (see sect. 3.3). 
Assuming the four-body modes to be dominant implies 



for each lepton flavor I. Moreover, the smallness of the total decay width implies that the 
decay vertex displacement may be detectable. 

3.2 Bounds from existing LHC searches 

The challenge of detecting stop decays for compressed spectra is all in the capability of 
reconstructing and identifying the soft decay products of the two stop decays. On the other 
hand, this is not the only experimental handle we have. 

The first problem is with triggering these events. The jets and leptons originating from 
the stop decay are too soft to be used to retain the events during on-line selection, given the 
(CPU and bandwidth) budget of the experiments. The only possibility is to detect these 
processes through the associate jet production {it* plus one or more jets), with a consequent 
reduction of the effective cross section. 

In the worst case scenario, all the decay products are lost and one is left with one or more 
jets: bounds exist from monojet searches [23, 24] (performed in the Dark Matter context) and 
from searches with > 2 jets [25, 27] (performed in supersymmetric contexts); see also [26]. 

In the best case scenario one can also detect the leptons from the decay of the stop pair. 
This is why it is interesting to consider a set of analyses that focus on > 2 jets for events 
with or without leptons. The CMS razor analysis [28, 29] is a all-in-one answer to our needs, 
with the additional advantage that the jet selection in the analysis is looser than the one 




(27) 
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Figure 7: Distribution of Mr (left), R 2 (center), and box-by-box event fraction (right) for pair- 
produced stop events as a function of the stop mass, for t — > tvgbN decays and m ( - — Mdm = 30 
GeV. Even if this case is the most favorable for the selection of leptonic final states, the hadronic 
box is the most populated due to the small value of — Mdm- 

used in the hadronic SUSY searches: pj? > 60 GeV for the first two jet; > 40 for the 
other jets. The looser jet selection increases the effective cross section we are sensitive to. 

To estimate the sensitivity of the search to the soft leptons from the stop decays, we 
implemented an emulation of the razor analysis, based on generator-level jets and leptons. 
We generate pair-produced stop squarks in y/s = 7 TeV pp collisions using PYTHIA8 [30]. 
The stop are forced to decay with a flat matrix element as t — > iugbN. The transverse 
momenta of all the visible particles are summed to compute the missing transverse energy at 
generator level. Similarly, these particles are clustered into jets using the FAST JET [31, 32] 
implementation of the antikT [33] algorithm. As for CMS, we use R = 0.5 to define the 
jet size. The razor variables and the six boxes (MuEle, MuMu, EleEle, Mu, Ele, and Had) 
are defined following the instructions provided by the CMS collaboration [34]. To take into 
account the limited efficiency in lepton detection, we applied the efficiency curves of the CMS 
dilepton SUSY search [35], using a hit-or-miss analysis. This is a valid procedure, since the 
lepton definition in the razor and dilepton SUSY searches are similar. 

We scan the value of the stop mass between 100 GeV and 400 GeV, fixing the stop-to- 
neutralino mass gap to 30 GeV. We show in fig. 7 the distribution of the razor variables for 
different stop masses, as well as the breakdown in boxes. A few important features should 
be noticed: 

i) The Mr variable approximates the momentum of the jets in the frame such that 
\Ph I = \Pj2 1 • I n the case of squark pair-production, for which this variable was designed, 
this corresponds to the squarks rest frame. This is why the Mr distribution for this 
case peaks at the Ma = (m~ — M^ M )/mf 1 . Instead, in the case we consider here the 
jets come from the associated (non-resonant) production and the peak is at ~ 150 
GeV, regardless of the stop and neutralino masses (due to the selection on the jet 
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and not to the SUSY kinematic). 

ii) The R variable is defined as M^/Mr where < Ma is a transverse invariant mass, 
such that the QCD background peaks at R ~ 0, while the signal can produce events 
with larger values of R, where the two jets have similar directions, opposite to the 
direction of the neutralinos. In the case of compressed stop spectra the R 2 distribution 
has some dependence on the stop mass, due to the correlation between the stop mass 
(setting the scale of the hard interaction), the spectrum of the associated jets, and the 
missing energy in the event. 

iii) The majority of the events selected by the analysis falls in the hadronic box. 

All these features are explained by the fact that the analysis is only sensitive to the events 
with two associated jets. These jets form the hemispheres and the razor variables are com- 
puted for a non-resonant production. In the largest fraction of the events the decay products 
of the stop are not seen and, effectively, the signal behaves like for the direct production of 
Dark Matter [36] , The stop plays the part of the Dark Matter, with the big advantage of the 
production cross section much larger than for Dark Matter direct detection. At the same 
time, the result is largely independent on the final state the stop decays to. 

These considerations suggest that the Had box is the only relevant sample to consider in 
our study. This is also the only box for which the information needed for phenomenology 
studies (observed yield and expected background vs R 2 and Mr) are given. While we limit 
our study to the Had box, we stress the fact that there is some sensitivity in the Mu box and 
Ele boxes, which could have be exploited if we had the relevant information. The importance 
of the Had box over the others also implies that the monojet analysis is a good candidate to 
look for our signal, as it is for Dark Matter direct production. 

We show in fig. 8 the limits obtained with the monojet and the razor (Had-only) anal- 
yses. For both the analyses, we consider the expected background yield (with error) and 
the observed yield, and we model the likelihood according to a Poisson distribution. The 
background uncertainty is described using a log-normal function. We assign a 30% error 
to the signal efficiency, to take into account the differences between our implementation 
of the analysis and a more realistic description of the CMS detector. We then derive a 
posterior-probability density function for the signal cross section as: 



where b (e) is the actual value for the background yield (the efficiency), b (e) is its expected 
value, and Sf, (8 t ) the associated error; Ln(x|m, 8) is a log-normal function for x with mean 




(28) 
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m and variance a; n is the observed yield, L is the available luminosity (for which we neglect 
the ~ 4% error) and a is the signal cross section. In the case of the razor analysis, the actual 
posterior is obtained as the product of the posteriors in each of the bins provided in [34]. 
We verified that taking L = e = 1 we can reproduce the limit on the signal strength for the 
two analyses. The 95%-probability limit is obtained integrating the posterior from up to 
the value ctup such that 

/ CTUP P(a)da 



0.95. (29) 



The left plot of fig. 8 shows the 95%-probability limit on the signal cross section as a 
function of the stop mass for both the analyses, fixing the mass split at 30 GeV and 100 GeV 
(with the stop decayed to t*N). The sensitivity of the monojet analysis is limited by the 
tight selection on jets and missing transverse energy. The limit is worse for larger splitting 
because of the veto on any third jet with > 30 GeV. At the contrary, the razor analysis 
is more efficient for this signature and more performant for larger splitting, since no veto is 
applied. One should also consider that at large values of the mass splitting the five leptonic 
boxes could further improve the sensitivity. 

The right plot of fig. 8 shows the limit in the stop mass vs neutralino mass plane. This 
plot shows the same qualitative features as the ID limit plot. At large splits, the limit from 
the razor analysis is found to be consistent with (and slightly worse than) the official limit on 
stop pair production [29]. Both the ID and 2D limits were obtained comparing the excluded 
cross section with the NLO+NLL tt* cross section at 7 TeV taking the decoupling limit for 
the other SUSY particles [37]. 



3.3 Dedicated analyses 

The existing limit is interesting, considering how challenging this signature is. This study 
also shows once more that the inclusive searches by ATLAS and CMS are much more general 
than the signal signatures they have been designed for. While a dedicated search could do 
better for a specific scenario, the inclusive searches are a good assurance policy for unexpected 
signatures. Repeating the analysis at 8 TeV with more data will certainly push the sensitivity 
further. On the other hand, we think it is interesting to imagine how the analyses could be 
changed to improve the sensitivity. 

One could certainly gain by using looser kinematic requirements. The limiting factor 
is related to the triggers. For example, it was pointed out extending the razor analysis at 
the tail of R 2 for low Mr could improve the sensitivity to DM production [36]. The same 
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Figure 8: Left: predicted cross section and experimental limits as functions of the lightest stop 
mass. Right: excluded regions in the (m^, MdmJ plane from ATLAS analyses (dotted regions 
shaded in yellow), CMS analyses (dot-dashed regions shaded in green), our re-analysis (red). 

conclusion applies to compressed stop-neutralino spectra, since the signature in the razor 
Had box is the same. 

A change in the lepton selection could further increase the sensitivity of these analyses. 
The left plot on Fig. 9 shows the distribution of the muon px for W^+jets events selected 
by the CMS monojet analysis, before applying the muon veto and the isolated track veto. 
This is compared to the equivalent distribution obtained for events with pair-produced stops, 
decaying to W*bN, with at least one of the two W* producing a \w pair. We consider two 
values of the stop mass (m t ~ = 150 GeV and m 4 ~ = 270 GeV) for AM = 15 GeV. Requiring 
one muon with px < 15 GeV corresponds to reducing the Z(uv) +jets background to a 
negligible level, and to rejecting ~ 92% of the other backgrounds. 

To evaluate the potential improvement due to this change, we applied the monojet anal- 
ysis to the generated stop-stop samples, and we separate the selected events in two boxes 
(as for the razor analysis): the Muon box, including all the events with one muon with 
Pt < 15 GeV; the Had box, with all the other events. We then distribute the background 
in the two boxes as follows: all the Z(vi/)+jets background to the Had box; 8% (92%) of 
the other background in the Mu (Had) box. We then evaluate the potential sensitivity of 
this modified analysis on a sample of pair-produced stop decays, decaying to W*bN, 20% of 
which produce at least one muon in stop decay. 



15 



= 0.7- 
n 

co : 
"§0.6- 

Q. 

0.5r 
0.4r 
0.3 r 
0.2r 
0.1 = 



W+jets 

1 1 +jets m_ = 250 GeV m LSp = 235 GeV 
1 1 +jets m. = 300 GeV m LSp = 235 GeV 




10 20 30 40 50 60 70 80 90 100 

muon p T [GeV] 



CMS Monojet Analysis AM = 15 GeV 



Modified Monojet+MU AnalysisAM = 1 5 GeV 




150 200 250 300 

m~ t [GeV] 



Figure 9: Improvements that can be obtained with a dedicated search. Left: distribution of the 
muon pt for W+jets and stop-pair events passing the CMS monojet selection criteria, except for 
the muon veto and the veto on isolated tracks. Right: expected excluded cross section for stop pair 
production obtained from the CMS monojet analysis (blue) and a modified monojet+muon search 
(black). Events are generated with four-body stop decays to ff'bN, of which ~ 20% produce one 
muon. 



The right plot of Fig. 9 shows the expected exclusion limit, compared to what is obtained 
with the usual monojet analysis. A similar improvement could be used for electrons, provided 
the understanding of the electron identification and the fake rate at low Pt- One should keep 
in mind that our results come from a simplified description of the CMS detector. A more 
accurate assessment of the improvement can only be obtained with a detailed simulation 
of the detector performances. We look forward to see this change applied to the monojet 
analyses by ATLAS and CMS. 

As a side remark, we would like to stress the fact that the stop decay products could 
be displaced from the primary vertex of the proton-proton collision. Requiring a displaced 
vertex, particularly with one muon originating from it, can potentially reduce the standard 
model background to a very small level. However for AM 30 GeV only a small fraction 
of the t\ decay after a detectable path of about 1 mm. An accurate estimate of the signal 
sensitivity for a diplaced- vertex analysis would require an accurate description of the vertex 
resolution for the LHC detectors and should be investigated directly by the ATLAS and CMS 
collaborations. Even if the signal reduction is too large to be beneficial for 8 TeV searches, 
this is an interesting possibility in light of the high statistics expected for the future 14 TeV 
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LHC run. 



4 Conclusions 

We have put forward a series of theoretical arguments that motivate the existence of a mainly 
right-handed stop in the mass range m~ tx = 200-400 GeV, together with a neutralino 30- 
40 GeV lighter and, possibly, a gluino with mass below 1.5 TeV. However, quite independently 
of any of these specific motivations, the search for stops nearly degenerate with the neutralino 
LSP is an important experimental task, necessary to cover possible corners of parameter 
space where supersymmetry may still hide. 

We have pointed out that, when the mass splitting between stop and neutralino is smaller 
than M w + m b , the previously-neglected four-body decay processes tx —> Nbt+vi and t\ — > 
Nbqq' can compete with the flavor-changing decay ti — > cN. The presence of a charged 
lepton in the final state of the four-body decay gives a useful handle to identify the stop in 
this experimentally difficult mass configuration where all visible particles are relatively soft. 

Regardless of the particular stop decay mode, the request of an extra jet in association 
with the stop pair greatly improves triggering capability and signal identification. We have 
shown that the inclusive searches using razor analysis are very efficient to probe stops nearly 
degenerate with neutralinos. In this region of mass parameters, we are able to set limits on 
the stop that are stronger than those published by ATLAS and CMS. Our limits (shown 
in fig. 8) extend up to stop masses of about 250 GeV, even for a vanishing stop-neutralino 
mass difference. This means that the LHC has already started to probe the "light stop 
window" motivated by our theoretical considerations, but most of the interesting region will 
be explored only at LHC 14. 

Note Added 

While our paper was being completed, similar results were presented in [38, 39]. In our analysis 
we used the response function for the CMS detector, provided by the CMS collaboration, instead 
of trying to emulate the LHC detectors. More importantly, we used the full likelihood provided 
by the CMS collaboration for the inclusive razor analysis, which gives a realistic description of the 
likelihood resulting from the data. This prevented us from extending our study to the razor btag 
search [40]. The latter has a better sensitivity in presence of bjets with pt > 40 GeV (i.e. far 
from the diagonal of the m^o vs. plane), if an accurate emulation of the btag efficiency and the 
mistag rate is reached. On the diagonal, only a small fraction of the associated jets are bjets, such 
that requiring a btag has the effect of reducing the expected signal much more than the factor-three 
background reduction. 
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Appendix 

In this appendix we derive the four-body stop decay width given in eq. (24). In the limit of small 
mass difference (AM <C Mw,mt), the amplitude for pure right-handed stop decay tR — > Nbe + v is 

KP = ^gf=(i* (30) 

where P{ are the quadri-momenta of the particles involved. The decay width is given by 

The 4-body phase space integral dcj)^ can be analytically performed at leading order in AM = 
JTTf — Mdm- Indeed, by writing the decay as t — > XY — > (Ne)(bis), the amplitude for each sub-decay 
is separately Lorentz invariant. Thus, using 

(4) = dsxd^^ d(f>{2) ^ ^^(2)^ ^ hv) (32) 

(2-kY 

2 8# 6 tan 2 6> w 2 

= QA / f4 TT)2 S x(sy-M DM ), (33) 



we get 



in- 



(mj-y/BY) 2 ^ M 2 M \ A(m?,s x ,s y )K| 2 2c7 6 tan 2 # w I 



r= ds Y dsx\l- *' U % = V , (34) 

A/2 M io V s ^ / 4(4vr)5m? 9(47r) 5 M^m 2 m t ~ v 7 

7=™}/', i */ <I -^ <b 5f,-^%y i /a+*-») > -*-w- (35 » 

yMg M /m? 7o y V mf / v 315 

where we have kept only the leading order in AM. From these expressions we obtain eq. (24). 
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